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SUMMARY 


This  report  documents  the  development  of  a  mathematical  model  which  describes 
the  gunner's  performance  In  an  AAA  tracer-directed  manual  firing  task  under 
periodic  observation  Interruptions.  Observations  are  interrupted  via  blanking 
the  target  aircraft  from  the  optical  display.  During  the  interruption 
period,  the  gunner's  performance  on  minimizing  tracer-to-target  error  is 
considerably  degraded.  Reduced-order  observer  theory  is  applied  to  design  a 
blanking  gunner  model.  The  model  consists  of  a  reduced-order  observer,  a 
linear  feedback  controller,  and  a  stochastic  remnant  element.  Both  the 
tracking  and  the  tracer  errors  are  considered  measurable.  The  effect  of 
observation  interruption  is  modeled  by  degrading  the  model  parameters  per¬ 
taining  to  the  observed  states.  An  exponential  decay  form  is  assumed  for 
these  parameters  during  the  blanking  and  recovery  periods.  Model  parameters 
and  associated  time  constants  are  Identified  systematically  from  empirical 
data  via  a  least-squares  minimization  algorithm.  Simulation  results  for 
model  predictions  versus  empirical  data,  over  several  blanking  conditions 
using  a  typical  helicopter  operational  trajectory,  are  included.  The  results 
show  that  the  gunner  model  can  adequately  describe  human  response  in  this 
compensatory  tracking  and  firing  task  under  observation  interruption.  The 
gunner  model  can  be  Incorporated  Into  existing  attrition  models  to  evaluate 
the  survivability  of  aircraft  in  tactical  engagement  scenarios  with  optical 
countermeasures  present . 
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PREFACE 


This  report  documents  a  study  performed  by  Systems  Research  Laboratories, 
Inc.  (SRL),  Dayton,  Ohio,  for  the  Air  Force  Aerospace  Medical  Research 
Laboratory  (AFAMRL) ,  Human  Engineering  Division,  Optical  Countermeasure 
program.  This  work  was  performed  under  Contract  F33615-79-C-0500.  The 
Contract  Monitor  was  Mr.  Donald  McKechnle  and  the  Program  Manager  was 
Maj.  Allan  M.  Dickson.  The  SRL  Project  Manager  was  Mr.  Kalle  Bishop. 

The  author  extends  his  appreciation  to  Dr.  Carroll  Day,  Messrs.  Walter 
Summers  and  Marls  Vlkmanls  of  the  Human  Engineering  Division  of  the  Air 
Force  Aerospace  Medical  Research  Laboratory,  Wrlght-Patterson  AFB,  for 
their  valuable  comments  and  discussions. 
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Section  I 
INTRODUCTION 


The  modeling  of  human  performance  In  an  antiaircraft  artillery  (AAA)  system 
has  been  extensively  studied  by  many  Investigators  In  the  past  decade  [e.g., 
Klelnman  and  Perkins  (1974),  Phatek  et  al.  (1976),  Kou  et  al,  (1978)].  Most 
of  these  works  dealt  with  the  modeling  of  human  response  In  a  simple  tracking 
task.  In  the  event  of  Interrupted  observations,  the  operator's  tracking 
performance  degrades  significantly  during  the  Interruption  period  and  poses 
an  appealing  modeling  problem.  The  author  tackled  this  problem  by  degrading 
•several  observer  and  controller  gains  In  the  model  and  proved  to  be  rather 
successful  (Yu  et  al.,  1980).  Efforts  were  then  directed  to  study  human 
response  In  a  manual  tracking  and  firing  task.  In  this  task,  the  operator 
(gunner)  directly  controls  the  gun  turret  and  fires  tracer  rounds  continu¬ 
ously  toward  the  target.  The  gunner  perceives  the  tracer  ending  position 
and  continuously  adjusts  weapon  pointing  In  azimuth  and  elevation  to 
minimize  the  tracer-to-target  error.  In  this  system  mode,  the  gunner  has  to 
play  both  the  role  of  a  tracker  and  a  lead  angle  computer.  The  conventional 
tracking  task  Is  greatly  complicated  by  the  Inclusion  of  lead  angle  estima¬ 
tion.  Wei  (1981)  developed  an  observer  gunner  model  which  treated  the 
tracer  Information  as  delayed  measurements.  The  Intent  of  this  paper  Is  to 
extend  the  author's  previous  work  to  consider  an  even  more  general  tracking 
and  firing  scenario,  l.e.  to  consider  a  tracking  and  firing  task  subject  to 
external  measurement  Interruptions. 

The  Interruptions  occur.  In  the  real  world,  through  various  electronic/ 
optical  countermeasures,  weather,  or  terrain  conditions.  In  this  study, 
extensive  manned-simulation  experiments  were  conducted  at  the  Air  Force 
Aerospace  Medical  Research  Laboratory  of  Wrlght-Patterson  AFB,  Ohio.  A 
typical  helicopter  operational  trajectory  was  used  In  the  experiment.  The 
trajectory  consists  of  three  phases.  During  the  first  phase,  the  target  Is 
standing  still  at  certain  altitude  and  Is  half  masked  by  some  terrain  con¬ 
figuration.  At  the  onset  of  the  second  phase,  the  target  pops  up  for  a  full 
unmask  flight  and  moves  horizontally.  Blanking  of  target  Is  administered  In 
this  phase  only.  Blanking  durations  range  from  1.5  sec,  3.0  sec,  8.0  sec, 


and  full  blanking.  Certain  repetition  of  blanking  duration  is  also 
Included. 

The  structure  of  the  gunner  model  in  a  tracer-directed  fire  system  In  Wei 
(1981)  is  adopted  and  generalized  here.  Nonlinear  ballistic  equation  Is 
used  to  compute  the  loci  of  elevation  projectiles.  The  model  consists  of  a 
reduced-order  observer,  a  linear  feedback  controller,  and  a  noise  remnant 
element.  The  remnant  function  which  lumps  all  of  the  random  effects  due  to 
measurement  noise  and  human  neuromotor  response  noise  Is  assumed  to  be 
Gaussian  with  Its  covariance  being  a  function  of  estimated  target  velocity 
and  acceleration. 

The  effect  of  observation  Interruption  is  modeled  by  exponentially  degrading 
the  observer  gain,  the  controller  gains  pertained  to  observed  states,  and 
the  bias  term  in  covariance  function.  Model  parameters  and  time  constants 
are  Identified  separately  with  respect  to  no-blanking  and  blanking  empirical 
data  via  a  least-squares  minimization  algorithm.  The  computer  simulation  of 
the  designed  model  shows  that  the  model  predicted  tracking  and  tracer  errors 
are  In  good  agreement  with  empirical  data  over  various  blanking  conditions. 


Section  II 

AAA  TRACER-DIRECTED  FIRE  SYSTEM 


In  a  tracer-directed  fire  mode,  the  gunner  perceives  both  the  tracking 
error  and  the  tracer  error  on  a  two-dimensional  visual  display.  The  tracking 
error  e^^  is  the  difference  between  the  target  angle  6^  and  the  barrel 
pointing  angle  9g.  It  is  also  referred  to  as  "lag  angle"  later  in  this 
report.  The  tracer  error  is  the  difference  between  the  target  angle  6^ 
and  the  projectile  ending  angle  9  .  In  the  simulation  experiment,  the  pro¬ 
jectile  flight  path  ended  at  the  range  of  the  target.  Each  tracer  round 
disappeared  at  this  point,  9^,  from  the  display.  The  detailed  description 
of  the  configuration  is  described  in  Wei  (1981).  We  will  briefly  summarize 
the  underlying  dynamic  system  in  this  report.  Figure  1  is  the  block  diagram 
of  an  AAA  tracer-directed  fire  system. 
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Figure  1.  Block  Diagram  of  an  AAA  Tracer-Directed  Fire  System 


At  any  given  time,  the  target  trajectory  input  9^  is  fed  into  a  visual 
display  device  and  combined  with  the  barrel  pointing  angle  9g,  as  well  as 
the  projectile  ending  angle  9  ,  to  form  error  signals  e^^  and  e^.  The  human 
operator  observes  these  error  signals  and  generates  a  control  output  u  via  a 
controller,  or  H-grip,  displacement.  The  control  signal  then  drives  the 
barrel  and  rate  control  plant  for  a  new  barrel  pointing  angle  0g.  Tracer 
round  is  fired  at  this  angle  and  passes  through  the  projectile  ballistics 
computation  to  obtain  the  projectile  ending  angle  9^.  The  task  of  the 
gunner  is  to  constantly  align  the  projectile  ending  angle  to  the  target 
angle,  l.e.  to  minimize  the  tracer  error  e2.  The  dynamics  for  the  elevation 
and  the  azimuth  firing  system  are  very  similar.  In  addition,  the  elevation 
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system  can  be  decoupled  from  the  azimuth  system.  However,  the  azimuth 
system  cannot  be  separated  from  the  elevation  system  due  to  a  coupling 
factor  cos  (0g)£j^  the  measurement  equation. 

By  introducing  a  state  vector  Xj(t)  =  [x^j^(t),  x^2(t),  x^^Ct)]^,  "T"  means 
"the  transpose  of,"  with  x^^(t)  =  ej^.j.(t)  - 

and  x^2  “  i  =  1,  2*,  the  following  system  and  measurement  equations 

which  represent  the  underlying  tracer-directed  fire  system  can  be  derived, 
see  Wei  (1981) . 

ii  ’  u^(t)  +  E^(t)  u^(t-T)  + 

and 

Xx(t)  =  C^(t)  x^(t)  1  =  1,2  (2) 

where 


*lf  not  otherwise  specified,  the  first  subscript  index  1  represents  the 
elevation  (1  -  1)  or  azimuth  axis  (i  =  2),  while  the  second  index  repre¬ 
sents  the  1-th  element  or  row  of  a  matrix. 


■  cos 

e^(t)  »  -1.34  X  (1-t)  X  [^1  +  (0.0052T+0. 0004861^)  sin  0iB(t-T)] 
e^Ct)  =  -1.28  X  (1-t) 

gj^Ct)  *  (0. 0052+0. 000972t)  x  t  x  cos  0j^B(t-T) 
g2(t)  =  0 


y^2  elevation  or  azimuth  components  of  the 

target  acceleration,  the  gunner’s  control  output  and  the  observed  tracking 
error  (lag  angle)  and  tracer  error,  respectively. 


If  we  Introduce  a  transformation  on  the  states  and  “  *^1^11* 

x^2»  *^3  “  then  Equations  (l)-(2)  can  be  rewritten  as  follows. 


^(t)  x^  +  B^(t)  u^(t)  +  E^(t)  u^  (t-x) 


+  F. 


0iT(t) 


+  G|(t) 


X.(t)  -  x^^t) 


( 


where 


and  2^(t)  =  that  the  coupling  factor 

cos  9j^g  Is  removed  from  the  measurement  equation  and  absorbed  into  the 
system  equation.  Equations  (3) -(4)  represent  a  nonhomogeneous  linear 
time-varying  system  with  a  time-varying  delay  in  the  control. 


Section  III 

AN  AAA  GUNNER  BLANKING  MODEL 


In  Wei  (1981),  the  author  proposed  an  observer  gunner  model  for  gunner 
performance  in  a  tracer-directed  fire  system.  The  function  of  the  gunner 
can  be  decomposed  into  two  parts  to  be  modeled.  In  the  first  part,  the 
gunner  observes  continuous  signals  and  makes  an  estimate  of  system  states 
based  on  his  Internal  model  of  target  motion.  In  the  second  part,  the 
gunner  utilizes  the  observed  and  estimated  states  to  form  and  exercise  a 
control  action  in  order  to  achieve  his  objective.  The  former  one  corre¬ 
sponds  to  an  estimation  process,  while  the  latter  corresponds  to  a  control 
process.  The  reduced-order  observer  is  used  in  conjunction  with  a  linear 
feedback  control  law  to  model  the  gunner's  function.  The  structure  of  the 
model  is  shown  in  Figure  2. 


Figure  2.  Block  Diagram  of  an  AAA  Gunner  Model 

This  model  structure  is  retained  for  the  blanking  case  except  the  ballistic 
equation  is  no  longer  parameterized  by  a  linear  equation  relating  and  0^. 
Instead,  a  more  realistic  nonlinear  ballistic  equation  is  used  as  shown  in 
Equations  (3)-(4).  In  addition,  time-varying  gains  are  used  to  model  the 
effect  of  observation  interruption.  We  will  discuss  the  no-blanking  case 
first,  then  the  blanking  case. 


MO  BLANKING 


The  equation  representing  the  gunner's  Internal  model  of  the  tracking  and 
firing  system  can  be  written  as 

i^(t)  -  A^(t)x|(t)  +  ^(t)u^(t)  +  Ej^(t)  u^(t-T)  +  G|(t)  (5) 

2^(t)  -  x|(t)  1  =  1,2  (6) 

Since  both  and  measurable,  the  only  state  that  needs  to  be 

estimated  Is  The  state  reconstructor  equation  for  derived 

by  applying  the  reduced-order  observer  theory  (Luenberger,  1971) 


,1^  >s  •  • 

-bikiiCiUi(t)  -  "  '^12'^!®!^'^^ 

The  objective  of  the  gunner  Is  to  minimize  the  tracer  error  so  that  a  maximum 
probability  of  hit  could  result.  In  other  words,  the  gunner's  response  in 
the  control  process  would  be  to  stabilize  the  underlying  system,  especially 
the  tracer  error  x^2^b);  therefore,  a  linear  feedback  control  law  of  the 
following  form  Is  designed  to  achieve  this  objective. 

UfCt)  -  x^  (t)  +  v^(t)  (8) 

where 


[y 


11’  '12' 


Yisl 


Is  a  vector  of  controller  gains  to  be  Identified, 


Is  a  vector  of  measurable  states  and  estimated  state,  v^(t)  is  a  remnant 
noise  function  assumed  to  be  Gaussian  with  zero  mean  and  a  covariance 
function 


[v^  (t)  v^  (s)]  -  l^iT  <*^>1  +  “i3  Kt^">| 


for  all  t  and  s.  are  nonnegative  model  parameters  to  be  determined. 

0^^  and  0^^  are  estimated  target  angle  rate  and  acceleration,  respectively. 
Equations  (7)  and  (8)  represent  the  gunner's  response  in  the  estimation  and 
control  process  of  the  tracking  and  firing  task.  If  we  define  a  new  state 
vector 

5^(t)  “  ~  *13 

then  the  state  equation  of  the  closed-loop  system  is  obtained  by  combining 
Equations  (7)  and  (8)  with  Equations  (3)  and  (4)  of  the  actual  tracking  and 
firing  system. 

^(t)  -  A^(t)^(t)  +  D^(t)  2^  (t-r)  + 

+  ^J^(t)v^(t-T)  +  R^(t)  (1( 


where 


■  -  jT  u-  .*  -  ‘V 


There  are  seven  model  parameters  in  total,  l.e.,  k^,  Yj^2’  "^13  “il’  “12' 

and  that  need  to  be  determined  from  the  empirical  no-blanking  tracking 
data. 

BLANKING 

The  optical  display  of  target  was  blanked  periodically  according  to  the 
duty  cycles  and  durations  listed  in  Table  1. 


TABLE  1.  BLANKING  CONDITIONS 


Condition 

Duty  Cycle 

a) 

Blanking  duration 
(sec) 

1 

25 

1.5 

2 

25 

3.0 

3 

25 

6.0 

4 

50 

1.5 

5 

50 

3.0 

6 

50 

6.0 

7 

75  * 

1.5 

8 

75 

3.0 

9 

75 

6.0 

10 

100 

1.5 
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The  duty  cycle  Is  defined  as  the  ratio  of  the  blanking  duration  to  the  cycle 
time.  The  blanking  duration  Is  the  length  of  time  that  the  target  Is 
blanked  so  that  the  subject  cannot  see  the  target.  The  blanking  always 
occurs  at  the  last  portion  of  a  cycle  and  may  reoccur  periodically  over  the 
entire  TOW  firing  period.  An  example  of  blanking  sequence  Is  given  In 
Figure  3. 

BLANKING 

PERIODS 


4.0  8.0  12.0  16.0  20.0  T(SFC) 

Figure  3.  Sequence  of  Blanking  for  Condition  4 

The  gunner's  performance  deteriorates  considerably  under  observation 
Interruption  via  blanking  the  target.  In  Yu  (1981),  the  effect  of  blanking 
on  the  gunner's  tracking  performance  was  modeled  successfully  by  degrading 
the  gimner's  estimation  gain  k(t)  and  the  controller  gain  Y(t).  A  similar 
approach  Is  adopted  here  to  model  the  effect  of  blanking  In  a  more  complex 
firing  task.  More  specifically,  the  observer  gain  k^  and  controller  gains 
and  y^2  which  pertain  to  the  observed  states  xJ^j^  and  assumed  to 

decrease  exponentially  as  the  blanking  starts  and  to  Increase  exponentially 
as  the  blanking  stops  (see  Figure  4) . 
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During  the  recovery  period:* 


ki(t)  -  + 

aii(t)  -  a^^(tp 

1  -  exp  (-  j 

1  / J 

1  - 


(15) 


(16) 


The  time  constants  associated  with  each  gain  parameter, 

determined  from  the  empirical  tracking  data  collected  in  the 
blanking  experiments,  as  shown  in  the  next  section. 


*In  the  simulation  program,  the  length  of  the  recovery  period  is  defined  as 
the  minimum  of  1.5  sec  and  one-third  of  blanking  period  to  avoid  covariance 
being  negative. 


Section  IV 

PARAMETER  IDENTIFICATION  AND  SIMULATION 


The  least-squares  Identification  program  developed  In  Wei  (1981)  was 
modified  to  Identify  the  no-blanklng  parameters.  Equation  (10)  can  be  first 
decoupled  and  then  approximated,  via  the  Average  Approximation  Method,  by 
the  following  ordinary  differential  equation  (Banks  and  Burns,  1978). 


V4(t)  -  N^(t)  W^Ct)  +  M^  n^(t) 


Xi3(t)  =  0j^^(t) 


Xi^(t)  -  -k^(t)c^  ^IT^*^^ 


(17) 

(18) 
(19) 


where 


biCiYi2(t) 


-1 


1 

0 

1 

T 

T 

0 

1 

0 

T 

’l  0 

■xii(t) 

0  1 

xi2(t) 

0  0 

Xj^2(t-T) 

0  0  _ 

x' 12(1-1) 

1 

T 
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(l  +  b^Yi3)  x^3(t)  -  b^c^Y^3X^4(t)  +  b^c^v^(t) 


% 


c^x^3(t) 


Cie^(t) 


Yl3Xi3(t-T)  - 


Yi3Xi^(t-T) 


Vi(t-T)| 


Cigi(t) 


IDENTIFICATION  OF  MODEL  PARAMETERS 


The  equation  which  governs  the  mean  of  states  Is  obtained  by  taking 
expectation  of  Equation  (17): 


Wi(t)  =  N^(t)  W^(t)  +  M^  i^(t) 


(20) 


where 


W^(t) 


(l+biYi3)CiXi3(t)  -  b^c^Yj_3X^^ 


^(t)  = 


CiX^3(t) 


Ciei(t)Yi3 


Xi^(t-T)| 


Cigi(t) 


The  first  and  second  component  of  represent  the  model  prediction  of 
ensembled  mean  of  tracking  and  tracer  error,  respectively.  On  the  other 
hand,  the  covariance  matrix  ^^(t)  satisfies  the  following  equation: 


P^(t)  -  N^(t)  P^(t)  +  P^(t)  Nj(t)  +  Lj^(t)  S^(t)  L^^(t) 

where 

E^(t)  -  E 


(21) 


21 


I 


biC^ 


0 

0 


cie^(t) 


0 

0 


a^l(t)  + 


- 

^“13 

e,,(t) 

0 

ail(t-T)  + 

+  a^3 

e  ,^(t-T) 


The  first  and  second  diagonal  element  represent 

the  square  of  the  model  prediction  of  standard  deviation  of  tracking  and 
tracer  error,  respectively.  Notice  that  time-varying  parameters  are  assumed 
for  kj,  y^2»  reflect  the  effect  of  blanking.  Since  the 

blanking  effect  to  the  estimation  of  target  velocity  is  predominantly 

expressed  through  degradation  of  k^,  there  is  no  need  to  consider  a  time- 
varying  Yj^3.  a^2’  “i3' 


The  steady-state  value  of  the  parameters  are  first  identified  via  a  least- 
squares  curve-fitting  identification  program.  The  reference  curves  to  be 
fitted  are  obtained  from  empirical  tracking  and  tracer  data  collected  in  the 
manned  simulation  experiments  without  observation  interruption.  These 
experiments  were  conducted  on  an  AAA  simulator  at  the  Air  Force  Aerospace 
Medical  Research  Laboratory.  Three  simulated  helicopter  trajectories 
ranged  1500  M,  2000  M,  and  2500  M  from  the  AAA  system  were  used  as  target 
trajectories.  Figure  5  shows  some  characteristics  for  the  1500  M  trajec¬ 
tory.  Let  and  *3^2 empirical  ensemble  means  of  tracking 

and  tracer  errors  and  Sj,j^(t)  and  corresponding  standard 

deviations.  These  empirical  means  and  standard  deviations  were  obtained  by 
averaging  and  computing  the  variance  of  the  empirical  data  from  40  simula¬ 
tion  runs  with  the  same  target  trajectory  and  the  same  subject. 


.-'i 
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The  parameters  were  Identified  by  minimizing  the  cost  function 

<v] 

defined  as  follows: 


1  =  1,2 


where  t^  is  the  initial  time  when  a  selected  tracer  round  reaches  the  range 
of  the  target,  t^  is  the  time  when  the  last  tracer  round  is  fired,  is  a 
positive  weighting  factor  chosen  to  be  one  in  the  identification  runs. 

The  direct  search  algorithm  developed  in  Wei  (1981)  was  modified  to  identify 
the  steady-state  values  of  the  parameters.  The  tracking  and  tracer  data  of 
the  helicopter  trajectory,  ranged  1500  M  from  the  AAA  simulator,  without 
blanking,  were  used  to  obtain  the  following  steady-state  parameter  values 
shown  in  Table  2. 

The  time  constants  associated  with  the  parameters  were  determined  empirically 
from  the  data  of  the  1500  M  trajectory  with  blanking  condition  5  and  listed 
in  Table  3. 
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TABLE  2.  STEADY-STATE  PARAMETER  VALUES 


^Insensitive  for  elevation  case,  no  degradation  Is  necessary. 


Notice  that  the  time  constants  for  t,  ,  t  ,  and  t  are  the  same  for  both 

k’  Yj 

elevation  and  azimuth  gunner  model.  This  Is  as  expected  because  the  gunner 
manipulates  the  H-grlp  for  elevation  and  azimuth  tracking  Indlscrlmlnantly 
with  respect  to  the  observation  Interruption.  On  the  other  hand, 
and  azimuth  case  are  considerably  greater  than  that  for  the 

elevation  case.  This  reflects  the  steeper  Increase  of  uncertainty  to  the 
target's  position  along  the  azimuth  axis,  because  the  azimuth  component  of 
target  acceleration  Is  much  higher  than  the  elevation  component. 

SIMULATION  RESULTS 


The  gunner  model  was  implemented  on  a  CDC  CYBER  175  computer  to  simulate  the 
man- In- the- loop  AAA  tracking  and  firing  task.  For  the  convenience  of  numeri¬ 
cal  computation.  Equations  (20)  and  (21)  are  discretized  Into  the  following 
form: 


ij,"  +  h"  r" 

^  -1  -1  M 


(23) 


.n  _n 


(24) 


where 


r  • 

I- 


•r 


exp[l!L(t  )  •  o]  do  • 

1  n 


•:  i  ' 


exp[N^(t^)  •  a]  io  • 


ii  -iiK  > 


■  2l<'n> 


A  -  0.06  seconds 


A  simulation  program  was  developed  which  uses  the  recursive  Equations  (23) 
and  (24)  to  simulate  a  closed-loop  AAA  tracking  and  firing  task.  Inputs  to 
the  simulation  program  are  the  time  history  of  range  and  acceleration  of  the 
target  aircraft,  the  initial  angular  position  and  velocity  of  the  target, 
the  number  of  blanking  intervals,  and  the  blanking  intervals  in  chronological 
order.  Outputs  of  the  simulation  program  are  model  predicted  mean  tracking 
error  and  its  standard  deviation. 

Simulation  results  are  shown  in  Figures  6  through  17  for  the  blanking 
conditions  3,  4,  5,  6,  7,  and  10.  The  solid  curves  in  these  figures  are  the 
empirical  data  which  are  obtained  by  averaging  the  results  of  40  experi¬ 
mental  runs.  The  dashed  curve  is  the  model  prediction  of  ensembled  mean 
and  standard  deviation. 


Figure  6  and  Figure  7  show  the  comparison  of  model  versus  empirical 
elevation  mean  and  standard  deviation,  azimuth  mean  and  standard  deviation 
of  both  tracking  errors  (lag) ,  and  tracer  errors  for  the  no-blanking 


case.  Figure  10  and  Figure  11  show  the  results  for  blanking  condition  5 
which  had  a  SO  percent,  3.0  second  blanking  occur  during  [11.01,  14.01] 
seconds . 

Of  particular  Interest  is  the  comparison  of  the  empirical  standard  deviation 
In  Figure  6  and  Figure  10.  The  effect  of  blanking  the  target  to  gunner's 
performance  Is  clearly  demonstrated  by  the  sharp  Increase  of  the  standard 
deviation  of  tracking  errors  during  the  blanking  period  [11.01,  14.01] 
seconds.  This  effect  Is  very  well  modeled  by  degrading  selected  model 
parameters  as  Indicated  in  the  model  prediction  curve  in  Figure  8.  Similar 
agreements  between  the  empirical  data  and  the  model  prediction  can  be  found 
In  Figures  8,  9,  and  12  through  17. 

These  figures  show  that  the  designed  gunner  model  can  provide  consistent 
prediction  of  the  gunner's  empirical  tracking  data  as  well  as  the  tracer 
error  data  for  both  no-blanklng  and  blanking  cases.  These  figures  also 
demonstrate  that,  for  a  given  AAA  weapon  system,  the  same  set  of  parameter 
values  and  associated  time  constants  can  be  used  to  predict  the  human 
tracking  and  tracer  errors  for  all  simulated  blanking  conditions. 

However,  due  to  the  fact  that  the  helicopter  trajectory  has  very  low 
elevation  axis  maneuvering,  these  parameters  may  only  hold  for  similar 
types  of  low  maneuvering  helicopter  trajectories.  Reidentification  of 
these  parameters  may  be  needed  for  other  high  maneuvering  trajectories. 

The  computer  execution  time  of  the  overall  simulation  for  an  18  second 
helicopter  trajectory  takes  about  5.60  cp  seconds  on  a  CDC  CYBER  175 
computer . 


Elevation — No  Blanking 


Figure  7a.  Mean  and  Standard  Deviation  of  Tracking  Error 
Azimuth — No  Blanking 


Azlmuth—No  Blanklpg 


Figure  8a.  Mean  and  Standard  Deviation  of  Tracking  Ern 
Elevation — 1.5  Seconds,  50  Percent  Blanking 


Figure  9b.  Mean  and  Standard  Deviation  of  Tracer  Error- 
Azimuth — 1.5  Seconds,  50  Percent  Blanking 


Figure  10b.  Mean  and  Standard  Deviation  of  Tracer  Error- 
Elevation — 3.0  Seconds,  50  Percent  Blanking 


Figure  11a.  Mean  and  Standard  Deviation  of  Tracking  Error 
Azimuth — 3.0  Seconds,  50  Percent  Blanking 
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Azimuth — 3.0  Seconds,  50  Percent  Blanking 
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Figure  I2a.  Mean  and  Standard  Deviation  of  Tracking  Error- 
Elevation— 6.0  Seconds,  50  Percent  Blanking 
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Figure  13b.  Mean  and  Standard  Deviation  of  Tracer  Error* 
Azimuth-*^. 0  Seconds,  50  Percent  Blanking 


Figure  15b.  Mean  and  Standard  Deviation  of  Tracer  Error 
Azimuth — 1.5  Seconds,  75  Percent  Blanking 


Figure  16a.  Mean  and  Standard  Deviati 
Elevation — 1.5  Seconds,  1 


Elevation — 1.5  Seconds,  100  Percent  Blanking 
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Section  5 
CONCLUSION 

This  report  summarizes  the  modeling  of  a  gunner's  performance  In  a  complex 
AAA  tracking  and  firing  task  under  pseudorandom  observation  Interruptions. 

The  highlight  of  the  task  Is  that  the  gunner  fires  tracer  rounds  without  the 
aid  of  radar  and  lead  angle  computer.  Furthermore,  the  gunner's  performance 
Is  greatly  hindered  by  observation  Interruptions  via  blanking  the  target  on 
an  optical  display.  A  blanking  model  Is  designed  which  consists  of  a  reduced- 
order  observer,  a  linear  feedback  controller,  and  a  remnant  element. 

The  gunner's  performance  Is  parameterized  by  the  controller  and  estimator 
gains,  in  addition  to  the  covariance  coefficient  of  the  remnant.  The  effect 
of  blanking  Is  modeled  by  degrading  these  gains  and  coefficients  as  a  func¬ 
tion  of  blanking  duration.  An  exponential  decay  form  Is  assumed  for  these 
parameters.  The  associated  time  constants  are  determined  from  empirical 
data  collected  In  the  blanking  experiment.  A  direct  search  method  is  used 
to  Identify  model  parameters  systematically  while  minimizing  the  least- 
squares  error  between  the  model  output  and  the  empirical  data. 

Computer  simulation  of  the  proposed  gunner  model  show  that  the  model 
predictions  are  in  good  agreement  with  empirical  data  for  various  blanking 
patterns  using  a  typical  helicopter  trajectory.  These  results  demonstrate 
that  the  model  can  adequately  describe  the  gunner's  tracking  and  firing 
characteristics  in  an  AAA  weapon  system  subject  to  observation  blanking. 

This  gunner  (blanking)  model  has  been  incorporated  Into  the  MTQ  series  of 
POOl/OBS  AAA  engagement  models  and  Is  designated  as  program  POOl/OBS  3/6B. 

This  composite  program  POOl/OBS  3/6B  can  be  used  in  the  evaluation  of  air¬ 
craft  survivability  and  performing  weapons  effectiveness  studies.  Documen¬ 
tation  of  POOl/OBS  3/6B  is  in  preparation  at  the  Air  Force  Aerospace  Medical 
Research  Laboratory  and  will  be  distributed  separately. 


APPENDIX  A 

LISTING  OF  PARAMETER  IDENTIFICATION 
PROGRAM  ELEVATION  CASE 


COMMENT.  *»OHL6EI.IOtID»Lr60295tCVlt  lOENT  EL  MT*  USE  HiLt  TR«J* 


PIN. 

AITACH|1*P-lt }HU6HELSUBJ33iCY«ltNR«l. 


011EN3ITN  ALFHA(n,PSICri  .  A  (71  >0(71 1  ALP  HANI  It 

(jt^Sf5N'3UMS(fcriTl717?^('7,7T;r  I7i7),xl  (7.7>.3{7.7)~.F|6),' ?(  6»r> 
CCiMON/ARRAV/XI  taoaif  S(iaOO)|ELOO(lboO)  tTXUODD)  tTSIlOOOI.RANIlOOD 


•  ■ n>  I  \ 

CCNNaN/S/C0.CEL(NSrPtNOtNiT0,IPT,V10,Y2a 


IMTEGER  QtSUElTtFR 


LbSD 

00  1  lat.7 


READ' .Rl.  HOI  7.TO,IPT.Y10O,Y10 

~«rNT"v?,-(er;  adin.to  ,  iPTTnff - 

FCRMAT  (IH1,“I.O.OF  PTS  •  -  .  I4t  2X  i“OROER«  *»I2. 

■Tr<<T“isiiT  n pe-  *,si2.5//'ix, “REAo^EVCTf-^rfr"^ oi m yi o ^  -,512 

C.5I 


0EL*9.(IT*1PT 


K»K1-XT 

kEtbci*^! 

c  i»t,Rn 


EZMNrOUMZt  AZSD.  OU  M3t' 


RLAOdi  *3t  T  ,CUM6«AZf  AZOt  AZa0«EltEL0t0UPltAZMN«DUM2i  AZSOiOJMS 


_TX(IT»«5UMi» 

iStlGIalUMS 

COSTlNUt 


55 


$6 


CCSTINOi 
coal.  34 

TTSTBTITM - 

KEAO*.  (ALPHA  (1)  t  Ial,7) 


OLOJaAJ 


I  :J ’» «;  I  m  «  u  !f 


FC^MATI-OALPHAa  G 12 .5.“ I“»618 .5."t“»6  12.5, 612. 5»“l- 


.9.  •  itoie.at  »  I 

PRINT  4itAj,  iT.sueir 

TO^HAn’’  Ja  ■.G12.5t^  irERATI0WSa~*rr5~t"'S(BmR«THWSa-», 1 5r 

00  100  Kal,7 


Mill  40^  J  t0 1 


CO  12  Lalt7 


r  luiiL  i:ii«Er  i«:j7r  f  m 


1F((U.G£. 2t.AN0.(L.LE.4n  GO  TO  12 

raTF=»i:i»HAF 

CCNTINUI 

CALL  IHTGULFHANrim - 

IF  (AJ.GT.OLCJ)  60  TO  20 


PRINT  40.ALPHAN 


I  If  n  ti!  I  lra.li  ni  9  n-ii  i :  1 4 1 


U«i«C«l 


■  :»c4f 


CC  i;  Nat, 7 


C  .NTlNUs 


m-IM 


6C  TO  25 


IF  IA(K).LE.1.5»  A(KI-0.0 


CC  33  Lal>7 


■I 


’  n^i  in*  I*  i‘  I 


CKal.O 


IF  (CK.  (E.O.  0160  TO  100 


00  32  Nal,7 

■SUlianm^XI  Iliinaa2- 
^UNEaSUlEaXI  (2,HI*a2 


XlaSQTTCSUHl) 


X3aXl/X2 


FOTNAT  (aija  a.GlE.S/*  ALPHRat*, 612. 5,4|*, 612.5t* ia,612.5* 

— ja|»;6ir.5,aiaT612;y.*t*«;i2.5ia-|a>G12.5/ -  - 

3*  XKlIa  •iGlE.S/a  XI(ll7XI(2|a  a, 612. 51 


f-  (UI«C.3 

F2(L,H)*0.0_ 

CC'inNUi 

cc  117  i»i,y 
00  iir  <»i,  7 
x;(  j«K>so.o 


OC  120 

sui3»sm5*(A6s  1 0  <  n_» » _ 

IF  (3JHT.L£.£FSI60  TO  tOOO 
CO  130  >4«lt7 


OC  110  J>lt7 
Z<Nt  J)«0(NI»  V<N«JI 


CC'UINUi 
OC  140  j>lt7 


OC  140  L»li7 
CC  140  <sjf7 


Cf^TlNUi 


^LIAsO  .0 
OC  150  JsltZ 


S0'14»S014*XI  «1,  Jl»*2 
51114330^3  150741 


LO  155  J»l,7 
VdtJI^XKlt  J)/SUM4 


KC'lNTt, 
e*<OUNT-l 
OiT  l7S  <«ltQ 
OC  160  L«lt7 


F2(a,7»3F(Cll»V1X,HI  ♦  F2<Q,m 
170  ccsTInuS 
Fn)»0.3 


cc  lie  1«1.7 


ntjuiSo 


200  SU15(0)«SUH5  (Q)  »8  (KOUNTi  Nl  4»  2 

- 51  »5URn  JURTlTin - 

CC  215  43lt7 


IF  IK3U1T.LE.7)  CO  TO  159 

^T^Tr♦l - 

5U9IT«0 


«• 

00  250  <«1>7 


.  1600  CALL  EXIT 
,  tK1 

SU990UTINC  IKTCCALFHAtEJ) 

. - C151ireir/5R1?Sf7SEHP  {TOO  OTTSEWnrit  OfflTEDOriO  0017TX 1 1000 ) ,  TS  (1 00  0»- 

,  l,1«U(1010)tThET(1000i 


ftf  c  C  fi  I  *1 '  R  VV  i  jn  L  i*n  M 'f9  4  J I '  Fi  #4' 1 


OliENSION  HC4ltP(4t  4)  f  PI  I  4(41  ,P2(4f  41  (  ALPHA!  71  (A  (4(41 

i  ( e<  i&r;  ?BriTi7CbitiT(  isiri  »Tir;i;i7r*inT  i«T4I  imiV(  cqc<  4i4)i0(4i~ 

2(X3ll0  01t  (X4  (10  00)(EOH(10  301  (EOOH(  10001 

ecvh/alTncFTa  (i(i)(Tfi»T7lEAil,ii,Eeiii»  ( (ZAitinidl  (tei'i  t  tiif 
UlTA  HT'l./ 


57 


c 

C  IMTI>L!ttTtON 


SaL«C9‘*2/0£L 


OC  1  jaltNO 


CCS(I>  Jt*0. 


F(t)aO. 


PUiliaO. 00  0  0256279 


I  u  n>r 


X3(  11*0. 000 


ECH(tlaX3(ll 


lF(«9G.'jT,-20a.)  StaEXPIARGI 


C  COMPUTE  4N0  STORE  STATES  X3  AND  X4 


9b  FCRNATibGlE.S) 


K2a«>l 

X  3C  XITOXS  rXKT»EOa  (XIO  *aEC~~ 
lF(AL>Ht(ll  .CO.O.i  GO  TO  b 


>  >1  j  #*911 0  mwM  ti 


GO  ''0  3 


C  COMPUTE  ANO  STORE  ESTIMATED  TARGET  VEtOCITV  AMI  ACCELERATION 


ECMIKK)aX3(KK)*Rb(KKI 


MlaND**2 


C  START  INTEGTATIOA  LOOP 


H  I*)ar29 

-ISET*0 - 

CC  190  KXaltNSTP 


0119  4  WA*i  4  4 


IFIKK.LE.MI  GO  TO  100 


IFdSET.NE.l)  GO  TO  101 

IPEATi  (TTir^XEMrTKKI )  CRm=TX  (HT  •»r - 

SSla(SQRT(P(  Ulll-SEMPdi  >o»2a< SORT (P (2 tZM -TST 1 II** 2 


iUi  »J.a01<l».2*rAU(K<)  IAUIKK)»*2)*ilNITHfeMKK-N>> 


ALlaALPNAIR) 

KWfOTI 

Kla«»l 

RAaNAA/TAOUfXl 

A2«l.>(TAU(Kll>TAU(KKn/0EL 

" 

Cfl^aCO'AL^ 

A(l.li>«C  OR 

A  <l,»)a-C0*ALPHA<3T 

A(E,NC1I«A(1,1|«AL1*A2 

A(?.NO)aA(l,  2)TJH»A2 

00  2  Ial.N02 

AUlfD’RA 

..  2 

A  (11.  Jlla-RA 

CONTINUE 

C 

c 

94 


CR1«C0*ALPHA  Ut*ALl*A2 
CR<»«bO^ALI>HA  (Al 
SCAL1*S:AL»(AL1*A2J**2 

T?rT2r2r«'XL*>Ha  (?r*scAxr 


IFCtHOOCKK,  100I.EQ.0).O%.  (KK.CQ.ll)  PRINT  99.  TtR  (llAYlOtSIE  AN, 
1  L4Rr(P(l,lii,SS0  "■ 

FCRMAT(;G12.!) 


C  COMPUTE  MEAnI  tracking  error 


DC  110  IBI.NCIH 
.  Qo  12(1  Jcl.NCiH 

.  »EA 

."ITD  ccntinOe  = 

>  no  CCHflNUE 

, - FTi»  a  <  r.'^tinrjF  xarniK » scireF  xinxKi - 

F(»t«K5(KKl*CR5R<X<»<KK-Ml-XJ<KK*Mn*fl.O  01*(1.>A?)P|S.2*0.9T2P 

. - C  T4UIKKn»CTSlTHrr(RK-Tm - 

.  00  130  lal.NCiN 

. - cc  iiiir  lai.z - 

..  0(II«C(I)«CA1NTII,JI*FIJ) 

.TVS — - 

.  HIIi«D(Ii 

« - fiIT)aO'. - 

.  130  continue 

,,■5  ^ 

c  compute  error  due  to  mean  tracking  error 

."C - - - 

SMSAN»SiEAN*  «N«H-XEMP(KK»)**2»<H<2I-TX  «P»»p*2 

..■e - 

..  C  COMPUTE  COVARIANCE  MATRIX 

-j - 

CGC(1,1»>  (At.FHA(5)»ALPHA(6)AABS<E0H(KKI  )*AL>4A(r)»AIStE00H(KK>n 

.. - i~»5txi; - 

IFKK.GT.M)  CQC<2,2)a(ALPHA(5i»ALPMAI6i»AB$(E0M(KK-N)) 

- i — ♦  Al RM A « 7 » •  AfS lEBOUTKK^^ T»  ) •  SC Al  1 - 

CALL  NULTIEAINT.CQC.NOIN.  Nl.Pl.lO) 

,  (.JLL  HULnETiP.IIDIH.HlTRZVTlfl 

.  OC  E20  Itl.NCIN 

. - c(r»?(r  j«i,'wnH - - 

P(ltJia>l(It  J)*P2<I,UI 


en 

^1)R]JTINE  HaT(EtPtLtUi«HirlR} 

- DTiniTBtrmTirnr.istisi  ,nai — 

CO  10  latiL 

- - im - 

.  00  10  K>ltt 

OC  5  J«liLltl. 

.  TPtP^TC'iP'JE  ( ji'Fnn 

.5  ii»ir»i 

, - infiT<=Tm>T - 

„  Hi<K)»T"MP 

^•ni - 6  KO  »T£WP - 

IFtMR.Ea.llRETURN 


TIC”?0  niTL 
00  20  <«ItL 


II»K 

270  IJ  JaitLltL 

TtiP»TE-1P*G(  J)*E<IH 

IS 

nsiT*L 

KK»«<-1>*LH 

L2=L-l 

, - rrTurrmrz - 

,  L3ar*l 

, - Dir~T0  - 

K1«II-1»»L*J 

. - ic<*rj-ii*m - 

.  Jii  H(I(1)  *H(K2I 

. - fTTJ - 

SC^POUTINE  OSCRT(NOIH«AtQEl.tEA,EAINT|NTl 

, - cnnsiiN  R(nieA{iiiE«rNrMi,ccEFTJo> — 

C  SETS  EA«EXP|A»OELI  «EAlNTaINT£GRAL  EA  0  T3  OEL 

. - UnHliTnTRFl - 

.  NMNall'NOXN 

. - immrrsi - 

ccepint»»i. 

. - DD“lT"I<riNT  PI - 

, - nrxcEmii*xrEv»coErtn*iiFFvoArtir - 

C  NT  MUST  Be  AT  LEAST  3 

„ - -  CAUL  3IAgfW0IHtEAIRTt~AtS0EFTllt  COEPIEn 

DO  60  La3,NT 

. - CALrTfotr{ATE»mrt-i«Jitf|ioiriCriii - 

IFd.EQ.NTiGC  TO  TO 

.  - 6B~CALl  DTAGnroTWTtAIHTTgrtrTlt^COEFCUl - 

70  OC  no  XIb1«N6«N01H1 

■. - r>iiir»;ArrnFiTa - 

..  BO  CCNTINUi 

- - 

SIB’DJTINE  DlAGINOIHf  AtStCltCZI 

„ - ontMSTOR-Am  f»rti - 

NCtMltM3IM*l 

- HK»ioTio»imni - 

.  NM«NDH-1 

- - 

IFICl.El.l.Oi  GO  TO  10 

, - nr  s- j-HTiHf  rom - 


KaJ^NNl 
OG  6  lai.K 
k  AdlaBlI) 

A  «III>A(II)»C2 

KCTU^N 

npj - 


•  2.i»e,2t-0.Q1796At-B«>2B511 

tVTiTriilt  .12 AAJJtTini SlB, .  2f3»«T'-f7iT^ 5E“3i  .  ir  3 02t-'S 


APPENDIX  B 

LISTING  OF  PARAMETER  IDENTIFICATION 
PROGRAM  AZIMUTH  CASE 


,  ireir.T  JOTI,  :'47 5  to  O  , — t760295THEIt  758-39  60 - 

,  HmP(ONI  . 

.  CJHHENT ;  cH-.rszic.TnarrBff^qjTcramoENr  *z  p«tr-usr-HEn-TRA;j* 

.  FT  N, 


iTTWEEHErSTJrJSTTCTrrfJr?*!"!; - 

AT  TACHtTAP£7,  THL6T<;ACcftSU3J33t:ral«  I0a(.76a295 

FRT3^AM  dot  (TAPEl.TAPEZ, input, OUPUT) 

CUFTSITN  ALFhAiriVPSKFI  ,A(7T70{7)  ,ALrRA1«7T - 

ulTrSilTN  SUf't.Tbl.E  (7  l,i/(7,7)  ,Z  (7,7)  ,XI  l7,7 »  ,  3  T  7 ,7 » ,  F  ( 6)  ,  • 

- tCTH’TK/AP'^AI  /-XM00TJTrsriOa0)7A2OO(10  0QT,EL3(rO0DJ7?ANtlOD 

i.TXdoaT) ,  rs  (10001 

CC1M3N/S/:0, LEL,NSTPiNDrNiT0rn’T7Y10,Y207AZ30 - 

lCTIIAl  cm, par 

INTFJFP  ■3,SUEIT,FR - 

L(^a 

- N  PA  Pi7 - 

NfA?lsN^AP-l 

LiC  1  lal,  NPAF  ■ 

CC  1  Jal,NPAi; 

“"KT,  J)  =  <r  (I^  JTaJT.O  ^ 

IF  (I.£T.  Ji  U(l,J>aXI  (I,  JIal.O 
- ZnTJTaT.  0 

bC,  JJa  j.: 

i~ — cc'JTiKur  -  - - - - - 

FRal 

- RETD',  XI,  naThTir,TPrrrtinnrro - 

PRINT  I,?,  Kl,  MjlN,Ta  ,IPT,Y10 

— FCTNyTTlHI,-R-0;or  PTS-a--,m,7XrORaE1t«  “,I?i - 

C  2<,"INrT  TI7£a  “,5 12 . 5// IX,  "RE  AO  Ei/ERY  -,I2>“  POINT-, -,Y1 

C.51  . .  .  —  - 

M»T0/I1.03 

..  - - 

DtL»T.3T»IPT 


2(8,7) 
0)  — 


■07“ 


Oa  “,512 


aETOd,  ,3)  (I  ,tUM4,A2,A70,  DUHl,  EL,  ELD,  EL  C0,3  JN2,  ELMN,  DUH3,i  LSD, 

- C  Ial,  <T)  -  -  - . . .  -  - 

liC  *  *  1*1,  K 

- RETCd,  ,3)  T.COML.AZ,  AZO,-DUMl-,TLittO,ELTJOV3JN2,'’LMN,  DUM3,EL3D 

IF{;jC(i,) 

-03 - FCRMAT  (12012.5) - - 

!,»  iF(*n;  (r-i,iFT)  .NE.  01  go  to  on 

Iha  d-D/IPT  *1  - - - - 

IF(  I  -t.N-:.  l>  GO  TO  ■•9 

AZTC*AZ) 

■09 - CONTI  NOE - 

A  TDdHtaOUMl 

-  -  X  (IH)  =3  JM2  - - —  - 

tLT( lM)a£L-ELMN 

- -  S  (IN)  *0  JNT  - - 

itA  KIHI  *7.5 

- IFtDONi,.lT;2877DrT-RAN(IH|-aOOT1«f7-<930-r»Tl-9^'JN!») - - 

.  Ma^AM(H)/jEL 

- IFtH.LI.N)  GO  TO  <»4  - 

iS"T*I3£T*l 

IFdSET.NE.l  )  GO  TO  44  -  - 

I*  a? 

■  - CONTINUE - - - 

R  (1?(  2,  -.6  )  (T  1ME,0UN5,TEL,  0UN7 ,7 ELSO, 1 0*  UKTI 

DC  4’  I  *1,K  - - - 

ret:  (2,  .6)  T  lMe,DU(15,TEL,  OUNZ.TELSO 

■  “  *  IF[CTc  ,*|  )  1,8,48  - 

at)  F  C»aU  ( -.5  12  .  5) 


63 


- - 

ktlO*t  OUPHA  (l)tI>liNPA«l 

_ _ _ 

AjaOaO 

IS  OC  8  IsliNPAfi 
1 

filllsO.] 

A(I)s2.~(i 
CO'ITINU; 
n»SUBlT«d 

y  28<AZTAU-(AZ0*V10BI 

“CM.TT'iTGUO'MJrriJT - 

CL1J«AJ 

psTSf  <*JTTJrL^HrrMpriflp»i,  mparj 

fC^MAT(-OALPMA«  G 12 .5, “ l“, G  1*.5. "I'f  6  12.5t“l~t  G12.  5.“l“i 
t  C12.5|-|“,Gi2.S.*i“tG12.5)  ' 

PRINT  41,  AJtlTtSUBir  _ 

nrt^iKsr f-" t«  etraticihs*  ■''Tr57*sT»iT-<rrTroTs«“-,T5i 

00  100  <«1,NPAR 
Ttjirr»S'j9'irTi 
OC  12  L’lfNPAR 

-ipTHAKTa«^L  fH**n  *ETi?i  - - 

lF(((..G-:.2I.ANa.(L.L£.4n  GO  TO  12 

-jrniPtin  nnr:n  to  ror-AiPif  a  wTa^rminn - 

12  CCHTINU- 

- ciCTTN  rc  (xtroa  r 

If  (AJ.n.OlCJ)  GO  TO  20 

- 0L1  J»  A 1 — - - - - 

PRINT  43, (AtPhANINPI •NP«1,NPARI 
PXniT  4i,  AJ,  I?,505TT 
Cl<l«S(<)  «£(K) 

- rRr*r**<Kr - 

OC  15  N«1,NPAR 

- jriPflrnii  •  ALP  nrimfl - 

i»  ccntinus 

- inriiOTCTrrTynTXTTTni - 

GC  TO  25 

75  rr<»«-.5»xrin 

IF  (tlK).LE.  l.Si  A<KI«0.0 

7s  CM1T75 

OC  33  L>1,NPAR 

'  T  t~i  riTT  -  r  ■  B.  Cl  en  Tii~yB - 


8 


..  40 


41 

11 


IT 


cmi.o 

■twrrfJuF 


IF  IGK.  4E.0.  OIGO  TO  100 


^Ui'l'TUIZsO.  0 


00  3?  N>1,NPAR 


Tc^naGmr?xni;m**7' 

SUIBsSUIZ^Xl  (2,M)*»2 

"GwriNai - 


TT 


33 


Xl«STRT(SUni) 

ssrsKm  Sim?  r 

X3SX1/X2 


TRIHT~3 1,  OLD  J,'l  ACPffAmPr.  WP-T;iirPAWrrXT73tT- 


FCIHAT  ou*  »,G12.5/*  ALPHA«(*,  GIZ.  5,*l»,  J12.5,»  ••iGlZ.  5 

T»I*7?l2V5,*i*761?.5,Pf*rJ12;5^,*l*,G12.5,*l»/ - 

3*  XT(l>a  *,012.S/»  XH1)/XH2)«  ‘,G12.5) _ 


GC  TO  110 
CQNTINU; 


bC  TO  11 
siijso.a 


■cc'iirmrwwnr 

OC  lis  i>ltNFAR 

TCtTiOTlI - 

F2(L,K>10.0 


ITS — irc^rrrjt' 

DO  117  JcltNFAR 

- CO"  1 1 7~<»  l.'NTin? - 

117  X1(J,K)30.0 

cc"i?(r“r*irNTjR - 

120  S(.133SU13»(Ae&(0III)> 


TF'CSJiliLE.U-blGO  TO~Iinir 
OC  110  '1*1>NPAR 
■xcno”!*  liliFAR 
ZIS.J|3:)(Ni*V(N«Ji 
c’^TI^iaz 

LO  IbO  KltNFAR 


IVO 


me  IbD  LUJNFAH - 

OC  IbO  <*JiNFAR 
~xiTJ7  tTilfi  rr.trrFTTTGrr 
cc^tinu  j 

""scHii^ir.a - 

LC  ISO  Jel.NFAR 


.  TSTi — STTfc*  smeFx  I  fi ,  jr»  »7 
,  SCIbsSaiT  (SUChI 

, - I.Cnsr“J*TiNFAR - 

.  ISS  «  (1|  JlsXKl,  J)/!>UH4 


IbS 


IbO 


TTTT 


■RC'JNT.-?- 
b*<OJ»)T-l 

mtrr75-<'«-iTxi - 

CO  ISO  LeltNFAR 

"FoiiFnj  ♦~xr(KDtjNT,xT*vn(7rr 
co'iriNJ? 

"DC  iro  'X^lTWFAIl - 

F20fM|»F(Q)*W<KfH»  ♦  F2(a.M) 

mcTTTNtn - - - 

F<U  =  0.0 

'ccvTiHur - 

DC  ISO  1<1»NFAR 

m  « 0  j NT ,  r )  ■  xTnccoNr,  1 1 -F'z  lOt  i) 

SUiS(0)30.0 

-VC~T3ff-'i*TfHFaR - 

SU-IStO)  sSUMS  I01*B«<0UNT,M»**-2 

- soistai^saRT  isuMstari - 

DO  21S  IHiNFAR 

- V  (OUNT.M)  »t  (KOJNTt  Mr/SOKTTttT" 

215  CCXTINU; 


mr 


2U0 


■KCiNTf'otnrrn - 

IF  (<DUir.Lc.NPAR»  GO  TO  159 

-17»IT*1 - 

bUTITrO 

FK«  1 - - - 

CO  2S0  XxIfNFAR 


■TKT*TI— 

OKt^C.3 

"AtKl*2vr 


1 


I|  (do)  .£B(16)|£61NT  (16)|  - «(  HI4)  .EAI  NT  <«i  t<>>(  '  <  <»>  •  CQCi  «>t  *i>  t  C  1 4) 
2,)iT(ipj)  ,X4  (10  00)  f  EON (10  00)  fEOQH( too 01 

~t'(i'Ti\/»c?Nct‘  (jr(i,i)”*8(i)) , (tA(i,i),£B  (nr;(“Ai)iT (itiitEBiNT (iii 
DATA  WT/1./ 


initialization 


CE*C03(-LG(KK)) 

AA“*-e£L»4LPHA‘(  1)»TB 
1F(A^3.i;T.-200.)  Sla£XP(AR6) 
x3(<l)*X3(KK)»tDO(<K)*D£L 
IF(AlaH4(1) .EQ.O.)  go  TO  4 

■Ti4(  irffis  i^TJrJTfRK  j  *€oo(Kicr*tsui - 

4C  T3  1 

x4rxir?it<r(KT?T»^cDfi<K)^Osi: 

CCNTIN'Ji 

t  EMI  <  Ha  )(TJ  O)  -  X  4CK  1 ) 

1  F( « .  G£ .  2)  E COH  (KK)  ■  (  EOH  (  KK) 'EO  N( KK -1 ) ) /OEL 

-hTrAT(X<)7DEt - 

IF(«.L£.M)  GO  TO  10 

TS£TirS£T*i - 

lFtIi-T.NE.lt  GO  TO  10 
DU)  i^r  t)'  CS 
h  (»)»r2)»C6 


P(»(»>a(0.0fl715e7*CB)*»2 

TP'A  MS  TM(  i  1 7  Cf-  XEFlP  dT) »  »^2>~ni(2  )2CB=T  XCID  "  2 

5S0»(S1Rr  (P(1(1))/CB-SENP(1))<'»'  2*  ( SORT  ( P(2,  2) ) /CB-TS  (  1)  )•»  2 

. .  ' 

coriNu: _ 


r6a-r>4E10*T«M£LG(KKI  I 
SC«L«(C3»CtfJ*‘2/0EL 

*T«T*orc - 

AZal.-iTtUtKlI-IAUiKKn/JCL 
■X1R*C0»»LPH*  121  *08 
Aa*tla-C0R»Tb 


A(*fM01l«*C0fi»AZ 

- »T*7NrjSAHt2)»A2 - 

A(2«>)aT9 

x - 

C  COHPUIE  TSAlSIfICN  MATRIX  EA  AND  ITS  INTEGRAL  EAINT 

X - 

CALL  0SaRT(NCIHtbt3ELtEBtEBINT.5l 

- cir%»cr(r*An»HA  tlt*cb - 

c 

X“STATlTIHTfGWTIOKXTKIP - 

C 


C  1  SlPTCPdill  l,SSO 

-99 - FTEMATTSGTE.X) - 

C 


CO  1?0  JaliNClN 

- OTiTinTi  AtAUt  JJ’urji - 

120  ccsriNUi 

TTO - CC'irtHUi - - 

F  aM(C3>CRL)*X3(KKi  aCRA*XA(KKI 

- FtrTAX3(KTrr*X8 - 

lF(«.Gr.HI  F(21aF(2l  ACR4*  <  XG  IKK'Ht -XST  KK-H> )  *  A? 

- 00'13(r  lainreiH - - - 

UC  140  jal,2 

- oniaori)  aeaint  ti,j»*Fu> - 

lAO  COXTINUI 

- mn  *  cm - — - - 

CdtaO. 

T30 - C  CUT!  NU£ - 

C 


C'CONPOTE'CCVAPTAKCE'TtATRIX - - - 

C 

- cc:  I ITI)  •  fACFHrC5T*At1»HArBn>-AB5fFt3H«KKmAL-PHA  tr  r*ASS  CEODH  IKK)  ) ) 

1  PSCAL 

- CCTt2t2l»ALP)m5)»5CAL*A2»»2 - 

IFIKK.GT.N)  COC(2,2|a  I  ALPHA  (51  aALPHAI6)«ABSIE0HI  KK>H)  >  *ALP  HA  171  • 

- 1  -3SrE33H(Klf-H)l)*SCAL*A2»*? - 

CALL  HULTiEAlNT.CQGtNOIHi  NltPlt  10) 

- CALr-PULT(EA,P,KOI?1,Nl,»2,lO) - 

OC  ??C  IdtNClN 

- CC  -gPO  JaliACIH - 

Pllf  JIaPKIt  JI*P2(ItJ) 

-ZZV - CC'ITtNUt  -  - - — - - 


)-TS  m’-IST)-)- 


X^'CaNPUTi  ERR3P  OOE  Tff  STANOARO-  OEVlATtON 
C 


TW-  CCITlNUt" - - 

£3>0£l«(SHEAI««MT*SS0) 

- return  '  . 

EM 
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OHEMSIlN  E<  1)  *F(1I  tG(16)  tH(l) 
ficTU  r«ifL 

ii»i  _  _  _ 

tC  lJ 
TtiP»0. 

LrT"J«I7t:i7T - 

1t1P»T£1P*t(J)*F(II> 

KK>(<>1)»L»1 
h(«>»T£MP 
10  GC«»*TiMP 

iFiM^.ei.irfrpfu^N 
GO  20  I^l.L 
OC?0  k*l,L 
Tt  1P»0. 

II»k' 

00  13  J>IiLliL 

- T  iP  *6 1  ^iTFETin - 

lb _ _ _ 

20  h«<K»«T£MP 

i.2aL-l . 

OC  31  l3lfL2 
(.3«r»l 

CO  30  j3L3<i. 
kl«(t-l»*L*J 
K2*» 

30  H«i»«H(k'2r 

tn 

- bXt»oirriNnraTrT<N3iH,AV0EL;EA,  ehni'.-nt  » - 

„  DIMENSION  A  (DtEAdltEAINTdl  ,0  CEF<30) 

X - $£TS  eA«tXPtA«oeLr,£ArNT«IKr£GRAL  E A~im“3Er 

„  NtINl«N  JIN*! 

.  NMNOH'NOIH 

.  NH1*NT-1 

.  tCEF'fNTt»l. 

OC  10  laltNTPl 
V  il«NT-I 

-  10  CCE- (II)=0£L*COeFdI»l)/FLOATa) 

,  •£ - NT'  MUSVEE"  AT~CEA3r“3 

u  CALL  0I4G(N0ln.LAINT,  A,;0EF<l)«COeF(2tl 

«  BE  6fl  l*37NT 

CALL  KUlT  (A.EAlNTtNDIHiNNtEA,  1) 

.  iFtL.iQ.NTioT  Ttrro 

..  bO  CA.L  CI1G<NUIN,EAINT,EA,1.0,COEF<LI) 
fll  00  P  dT  ll»l  VNKf  NDlMl 
EA(II)>'AaXI>1.0 
fJT  CT^TTinJE 
..  EKO 

„  SLlffOWlNE  OlAGiNOITif  ATBt  CitC2} 

OIMENSION  A(l)«Bd) 

NLlPliNOlnVl 

NMNOIH'NOIM 

. - NRT=T;TMiT: - 

.  1I>1 

.  lF{Cf7E17f7ffJ  "GO  7a~ni 

.  _ LC  5  J*l»_NNi^mj1 _ 

.  ksJ^^Hl 

,  OC  A  IsJtK 

. - ir~A  (iTin  ‘sTn - 

.  A(TI)*A(IIi»C2 

.  ^  H«II*N)IM1 

. _ fl£’’U^N _ _ 

10  DC  r'j«l,NN,TNOlSi 
Ka J»Nni 
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APPEHDIX  C 

LISTING  OF  AM  AAA  GUNNER 
MODEL  SIMULATION  PROGRAM 


M8M,T20>CH70000.  1760 ! 95 tMEIt 2iS 3960 
CONHENT. •NEHOHLST'IUSi I O'L 760295.  :v«l* 

CO MHENT.  •AAX  HOOE6  3Law<I'IG  SIIJUATION  P90SB8M* 

AT  t'ACH.TAPE1.0HI.6HEL5!IBJ33tIO«L7  50  295.  CY«1  .MRS  1. 

FTN. 

(.30. 

PROGRAM  S  IHJGaMPJT.  OUTPJT.  T  APEl » 

COMHON/S/CK2t  .'1EL.IH.MOtM.riO(2i.X.T(2I.EL.ELO'l.  A200.KTAU.RA.  A2.NO 

_ Al.ND2.NAA.U;L.UAZ,ELTR.>2rT.lSEr  .21.22.  TAU.TSUSt,  T  £115)  .T.IBL 

C  ' 

C  THE  PURPOSE  OF  THIS  PROGRAM  IS  TO  SIMULATE  AN  ELEVTN  AND  A2 THUTH  TRA^ K 
C  TASK  IN  THE  TRA  lEP-O  IPECTEO  FIRE  (MODE  61  SYSTEM 
C  SUBJECT  TO  OPTICAL  3LAN<1N5 

C  INPUT!  THE  ELEVATtCN  (EL)  t  AZIMUTH  (A2)  ANGULAR  ACCELERATION  OF 

C  TARGET.  AND  BLANKINC  DURATIONSiUP  TO  15)  IN  CHRONOGICAL _ 

C  OROER 

C  OUTPUT!  MEAN  ANO  STAND  0£W  3^  LAG  AYGLE  _ 

C  •*•••  ALL  ANGLES  ARC  IN  UNITS  3F  RADIAN  *•••• 

C  TAUl  DELAY  IN  SECONDS 
C  ALPHA!  PARAMETER  VECT3R 

C  ELERR!  MEAN  EL  LAG  ANGLE  (I.E.  TARGET  ANGLE-BARREL  ANGLE) 

C  AZERRI  MEAN  AZ  LAG  ANGLE 

C  ELSDI  STANOARO  DEVIATION  Of  ELEVATION  LAG  AN6L|_ 

C  AZSOI  STANOARO  DEVIATION  OF  AZ  LAG  ANGLE 

C  ELTR!  MEAN  EL  TRACER  ERROR  (TARGET  ANGLE-TRACER  ENDING  ANGLE) 

C  AZTRI  MEAN  AZ  TRACER  ERROR 

c  elbar!  mean  el  rarrel  angle _ 

C  DEL  I  TIME  STEP  USED  IN  THE  INTEGRATION  ROUTINE 
C  TS(I)!  STARTING  TINE  OF  I-TH  3LANKIM6  DURATION 
C  TE(I)!  ENDING  TIME  3F  I-TH  BLANKING  DURATION 
C  YlOdl!  INITIAL  GUESS  OF  EL  LAG  ANG.E  _  _ 

C  Y102I  INITIAL  GUESS  OF  AZ  LAG  ANGLE 

C  UEL!  EL  CONTROL _ 

C  UAZI  AZ  CONTROL 

C  COTDI  EL  RATE  CONTROL  COEFF  _  _ 

C  C0(2)!  AZ  RATE  CONTROL  COEFF 

C  Kll  MO  OF  POINTS  IN  TNS  ENTIRE  TRAJECTORY  _ 

C  K!  NO  OF  POINTS  AFTER  T4E  fjrsT  TRACER  ROUND  IS  FIRED 

C  ELOD!  EL  ANGULAR  ACCELERATION  3F  TARGET _ 

C  AZOOl  AZ  ANGULAR  ACCELERATION  3F  TARGET 

C  X3I1)1  EL  ANGULAR  VELOCITY  OF  TARGET  _ 

C  X3(2)l  AZ  ANGULAR  VELOCITY  OF  TARGET 

C  x»!  ESTIMATION  ERROR  OF  ANGULAR  VELOCITY  OF  TARGET _ _ _ _ 

C  EL!  CL  ANGULAR  fOSITION  OF  TARGET 

C  M(l)l  MODEL  PREOICTEO  LAG  ANGLE _ 

C  M(2')i  MODEL  PREOICTEO  TRACER  ERROR 
C  Pd. II!  VARIANCE  OF  PPEOICTEO  LAI  ANGLE 
C  P(Z.2)I  VARIANCE  QF  PREOICTEO  TRAGER  ERROR 
C  TO!  THE  INITIAL  FIRING  TIME 
C 

_ REAQF.Kl.TO.IPT.  IPL _ 

PRINT  3.K1.T0.I»T 

3  FORPATdMl.-NO  OF  PTS»  ".  I  A.  2X.-INIT  TIME-  "  ,G12.5//1X.“REA0  EVERY 
C“.I2."  POINT"/) 

IFdBL.GT.OI  READ*.  (TSIO  .  TEI X)  .  K>1.  19L) 

PRINT  11.I9L 

11 _ format  dXj  I5j  IX j  "FLANKING  INTERVALS  ARE  "/i 

IFdBL.CE.i)  PRINT  i..lTS(<I.TF|i(|,K«l.IBL) 

I.  FORMAT  (5(1X.“(".  F9. 2.".  ",  F  J,  2.  ~)  “) ) 

KT«T0/0.03 
K»K1-KT 
T«T0 
MOlMaA 
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IPRINT*20/t!»T 

C3(l)<1.3<i 

C0(2>=1.29 

NDl>NOIM>t 

N02»N0IM-2 

NA»sN0IM/2-l _ 

rsrsi 

DEL=0.03»I“T 
CLSO*.OOS  "I.? 

A2SO».005»*1.5 

Z1=0. 

_zi»iu _ 

1SET*0 

UEL>0. 

UAZsQ. 

PRINT  7 

F0RHAT(/1H  ,2X,-TIMc“,9)(,“EL  XE.  - 1 9«  ,  “EL  ERR  5  X,  “E  LSD” ,  6Xt 

1“EL  CTR".6 X."A2  tf EL“t6X. "A ZERR“i 6X. “A2S0“.6X .“A2  CTR“ _ 

2t6X,~EL  TR",6X,"A?  TR"/» 

READ  (li  21  «T,OU^l|AZiAZa,AZOD«EL«ELO.ELOO«AI*(N,X,ZSO,S,I<ltKT 
C  > 

00  S  IsltK 

REAO<1.2)T,3UMt,AZ,AZD.tZl],ELt'LO.iLOO,AZt1N,X,ZSD,S 

IF(EOFtl)>lil _ 

IF<HdDll-l,IPT).NE.3)  33  13  5 
IHa*!-!*/!®!*! 

T*TO*<rH-H»OEL 

TAUsT.S 

IF(0UH1,LE.28TT. I  TAU*QUN1/ < 930. - . 19*00« II 

I1TAU»TAU/0EL _ 

IFCIST.EQ.tl  OTAUiTAU 

IST*IST*1  _  _  _ 

RAsNAA/TAU 
A2»1,-(TAU-0TALII /DEL 
OTAU*TAU 

X3<H-EL0 _ 

X3(2l«AZ0 

IF(<IM-1I  .Si.NrAUI  G3  13  9  _ 

TAUR»TAU 

IFCOUHl.LE. 41.011  T  AUR>r  AU*  ANAXK  0 . 6,  OUNl/5d3  O.i  _ 

nO(l)  «-TAUR*ELl-.  001*«5.2»TAUR*.4e5*TAOR**2l*COS<EL*b.’0  5) 

VIO 2—0.0  25 »SI3>H1..AZ0> _ 

Y102*-TAUR»AZ3 

IF<IM.N£.l)  GO  TO  13  _  _ _ _ _ 

ELO«EL 

AZO»AZ  _  _ _ 

Eio*yiO(i) 

-iiaami _ 

GO  TO  10 

ISET»ISETn  _ 

IFdSET.NE.ll  33  TO  10 

Zl»EL-(EL0-;10)*.101*«5.2*rAU».486*TAU»»2)*COS(EL0*0i05l _ 

Z2:AZ-IAZ0-E20I 

CALI  09SEt,6(EL£'^R.ELSDI _ 

ELBAR«EL-EL€RR 

IFIIIM-1) .Li.MTAUI  yiOC2)«ri02*:OSIiLeAR) 

CALL  OaSAZSIATERR, AZSO.EL3AR) 

IF((MOO<rH-t,I»?lyT),E3.0l  .3R.IIH.e3.1>»  PRINT  6,T,  X3«ll_»  ELERR, 
l«CLSD.U£LtXl(2I.AZERRtAZS0,UAZ.ELTR,AZTR 

FORWAT<UC12.SI _ 

C0NTINU£ 

FORMAT (12G12.5L  _  _  _ 

STOP 

ENO 


SU9ROIJTINC  39SEL6(ELERRtELS0) 

COHnON/S/CO(2t .ICL tKK.SO.y 10(2I«  X30I  21 .EL. EDO. A  ZOO .H.RA.AZ^ 


aoo 


A(N01tN02tN4«.U,UA?.;LT^ttZrR,IS£T.Zl tZ2.T£U, TS(IS), TEIISI .T,IBL 
OINEHSION  H(li),ACfc|l>>«RlUtb)»P2(i..l«>,  ALPHA(r).A(lt,  A) 

tE8{lS)tE3I^T{16)>£AiA,hl>EAlAT(A,At.F<A>  .OC  (A,A> ,  C(4l 
2«X3(1000)  «XAI1330)  .cDHdOOll  , £D]H(  10 OCI  . THET  (13001  •ALPSOCIOOO) 
EQUIVALENCE  I A  ( 1 , 1 )  ,3(11) .  (•  A  t  It  II  t  :B(  1)  ).  (E  AIHTd ,  1»  .EMKT  ( 1) ) 

DATA  MT/1./ _ 

DATA  ALPHA/1.  SAT  It  .  017491 1 . 024A33t  .42310.  .224A6E-7«  *179r5e-5t 
A  .17302E-3/ _ ; _ 

INITIALIZATION.  _  _ _ 

IF((KK-1).G;.H>  S3  TO  S _ 

IF(KK.GT.1I  SO  TO  5 

N1»N0**2  _ _ _  __ 

NOIHaNO 

ALPl>ALP10aALP-ll  111 _ _ _ _ _  _ 

ALP2>ALPHA(2) 

ALP3tiALPHAC3) _ 

ALP4aALPHA(4) 

ALPS>ALPHA(SI  _  _  _ 

ALP6*ALPHA(SI  '■ 

ALP7aALPHM7) _ _  _ _ 

ALP5D(1}>ALPS 

ACC«0. _ 

IFLAGaO 

SCAL«ca<ll**2/0EL  _  _  _ 

DO  1  I-ltNO 

00  i  J»1,N0 _ _ 

PdtUlaO. 

CQCd.JMO. _ 

I  AlItJIaO. 

DO  11  1*1, H9  _  _ _ _  _  _ 

N(II«0. 

O(II«0.  _ _  _ 

F(I)«0. 

II  PdtHaO.OOTO _ 

P(1.1I>0.00102S6279 

P(2, 21-0. 00003396711 _  _ _ _ 

X3lll>X3a(ll 

X4(l)>0.  _ _ _ _ 

EOH(l}>X3(ll 

_ EDDH(1I«0. _ 

S1>0. 

5  MdlaVlOdl  _ _ _  _ 

6  IFdSET.EQ.l)  4C2I«Z1 

IFdBL.LT.tl  53  TO  tO_  _  _  _ 

ALP5«ALPHA(5) 

1S»1FLAG*1 _ 

00  12  I*IS. I3L 

IFCT.GE.TSdl.AVO.T.LT.r-dl)  G3  T3  IS 
.TT>ANINtd.StA:C/3.> 

rFtT.GE.TEdl.ANO.T.LT.ITEddATTM  GO  TO  16 
ZFd. GE.dEddATTM  GO  TO  21 

_ ACCOj _ 

IFIT.LT.TSdSII  GO  TO  19 
GO  TO  12 

L6  ALPl«ALPlC*(AL»HAd)-AL»lJI  *d.-EXP<-C.4  3*(T-TEd)  III 

ALPS»-0. 0001*  d. -EXP  (-0.43*  (T-TS  dim 
60  TO  IS 

1 5 _ AC^  ACC»DEL _ 

xKag«i-i 

ALPl»ALPHAdl*S<P(-0.075»»(T-rS(ll  II 
ALP3»ALPHAdl*E<PI-9.52*(T.TSdl  n 
ALPS«0.0001»d.-EXP(-9.12MT-rS(im  I 
ALP10*ALP1 
W  ID  46 . 


21 


ALPt>Al.PHA(t) 

ALP3>ALP-|A(1I 
ALPS*ALPHA<9> 

CO  TO  IS 
12  CONTINUE 

18  CONTINUE _ _ 

ALPSO(KK)«tLP5  ' 

ARG»-DEL*AL'»1 

IF(ARG.GT.*200.»  S1>EXP(AR:i 
THETIKOsEL 
K1«KK*1 

K2»i<K-l _ 

C 

C  COMPUTE  TARGET  VELOIITV  ANO  ESTIMATION  ERROR 
C 

XSCKDsXS*^)  ♦i'JO’OSL 
ZFIALPl.EQ.a.)  so  TO  4 

_ X4(K1)«S1»X4«K'<I  ♦gOa»<l.-51>/AL»l _ 

GO  TO  3 

4  X4<K1)»X4<X<»*ET0*0£L  _  _ 

3  CONTINUE 

C0HIKK)aX3(<KI-X4l«)  _  _ 

IFCNK.GE.ZI  £0’)’t<KKt*<E0H(  <K) 'EOHi KZ » I /DEL 

_ X30(l)«X31Xlt _ 

C0R>C0(1)*ALP2 
CR4aC0(ll  •ALPA 
IFdKK-l)  .LE.NI  SO  TO  150 

ALlal.+0.(l0l*«5.2*TAU*0.481?TAy»*2>*SIM<_THET  ) 

All.lla.cOR 

A<li2l»-C0«t>»ALP3 _ 

A(2«N01)aA(l,lt<ALl«A2 
A<2fNO)*A<lf2)»ALl*A2_ 

00  2  Isl,N02 

J1»I*2  _ _  .  _ 

A<Jl(IlaRA 

A<J1.J1)««RA _ 

2  CONTINUE 

CALL  0SCRT«N0«S,0EL|EB,|B1XI,5)  _  _  _ 

CR3*CR4*AL1*A2 

SCAL1»SCAL*<AL1»A2»**2  _  _ 

C 

C  COMPUTE  MEAN  TRACKING  ERROR  tI.E.  LAG  ANGLE! _ 

C 

U»ALP2*Ha)  ♦AL»T*X<ZI»ALP4»(X3«KK1^<4JJ<K»J _ 

00  110  I«1.N0 

00  120  J»lfNO  _  _ _ _ _ 

0(I>>0(II«^EA(I,  J)*NIJ) 

120  _ CONTINUE _ 

110  CONTINUE 

F(l|a<l.-CR4)  *<3<KKI*CR4*X4IKX) 

F(2)«X3IKKI  ♦CR3*  CX4(ICK-'1I  -<3 « KK-MM  fO.  00  1*<  1.-A2J  •  »  5,2a  0.  972*1  AU 
^‘COSJTHETCCK-H  )  __  _ 

00  130  I>1,N9 

DO  140  J»l,2 _ 

0«l)iO«I)*EAlNT(I,J»*FU» 

140  CONTINUE  _  _  ^  _  _ 

MII)-0(I> 

OlllaO.  _  _  _ 

130  CONTINUE 

C _ 

C  COMPUTE  COVARIANCE  MATRIX 

C  _ 

CQCd,  ll>(ALP5*AL>6*AB5(E}-l(  Kk)l  AALP7*ABS(E00HIKK)  )) 

1  *50X1 

COC  C2,2»»(ALP«1C  K<-'1)*AL»S*A3S(£0M»<K-H»  »*ALPr*ABS«EOOHtKX-M) 

It  t  *80441 


call  MULrt^AiNr.c^CtNn, '41. 31,101 
CALL  HtlLT(*A,P,<IO.Nt,P2.1]| 

DO  22u  I«1.SD 

00  220  Jst.-^O  . 

PII,  JMP1(I.JI»32II,  J> 

22  0_  CONTIHUE _ 

150  CONTINUE 

ELERR=NI1»  __ 

ELSO^SORT  |P(1,1) ) 

ELTR«H(2) 

RETURN 

_ ENO _ 

SU3R0UTINE  3BSA76( AZERR.AZ30,£L:I 

CONNON/S/CO(2i  ,}EL«'«.NO,riO(2),  X30(2i,EL.ELOa.AZOO.H,RA, 

1  AZ.N01.N02.NAl,UEL.U,ELTR,AZrR,  ISET.Z1.Z2.TAU.  TS(  1 5)  ,TE  ( 15)  .T.IBL 

OIHENSION  NC4),3(l,,it),oiC4,4l,P2l4,itl,ALPHA(7>,  A(4,4)  _ 

1.8(16)  ,EB(15)  ,E5tNT(ie),£4(4,4l,EAINT<4,4),F<4l  .COC  (  4,4 )  .  0(  4) 

_  2.X3(10  00)  .X4(U0  0I  .  EOHtl  001 )  .  ED3H(  ID  00)  .  ALPSOtlDOC) _ 

EQUIVALENCE  ( A  ( 1 , 1 )  ,  3(  1) )  ,  (  EA  ( 1,  1)  ,  •  8(1)  ).  (EAINT(  1.  1)  ,EBI  KT  ( 1) ) 

DATA  HT/1,/  _ 

DATA  ALPNA/6.5T<)4,,1T«94..  1ZZZ3, 1.03  53,.  252S6E>5..  19766E>  3,.  757  651 

A-3/  _  _  _ _ 

IF((KK-1)  .GE.M)  GO  TO  6 

IF(<K.GT.l)  SO  ro  S _ 

NlaN0»«2 

C 

C  INITIALIZATION 
C 

00  1  lal.NO 

_  _  00  i_JaliNO _ 

P(i, J)«0. 

CQC(I,J)aO.  _ 

I  A(I,J)»0. 

00  11  1*1, NT  _  _ 

N(I)»0, 

0(I)*0. 

F(i)«0. 

II  P(I,  11*0.0100 
P(l,l)*(0.'lJ5918  3*COS(iL3)  I  **2 
P(2,2)»(C.n71517»C0S(EL;))**2 
X3(l)*X3a(Z) 

X4(l)»0. _ 

EbH(l)*X3(l) 

E00H(1)*3.  _ 

S1*0. 

ALPlsALP10*ALB'(A(l)  _  _ 

ALP2*ALPHA(2) 

ALP3*ALPHA(  () _ 

ALP4*ALPHA(4) 

ALP5>ALPHA(5) 

ALP6*ALPHA(5) 

ALP7»ALPNA(7) 

ALPSOdlaALOO 

_ ACC«0. _ 

iFLA6«0 

C  _ 

C  CONFUTE  ANO  STORE  STATES  X3  AVI  X4 

C  _  _ 

5  N(l)>riO(2) 

6  _ IF(ISET.EQ.l)  i((2)«Z2*;3S(iL6) _ 

CONfiNUE 
C9aC0S(ELG) 
tF(IBL.LT.l)  G3  TO  1* 

ALP5*ALPHA(5)  _ 

IS«IF».A6*1 
09. 12  t*IS,tBL 
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to  O  O  U  ^  OOO  ooo 


IF(T.i;t.TSII)  .AMO.!  .UT.rE(  in  33  T3  15 
ATT*AHINl(1.5,4:C/3.) 

IF(T.Ge.TE(U.AN0.T.LT.(r5(ri*ArTn  GO  TO  16 

IFtT.G£.(TElH*iTTH  SO  TO  21 

ACC^O. 

_ IFIT.LT.TSlISn  GO  TO  18 _ 

GO  TO  12 

15  ALPl»ALP10♦(AL'•^A^tl-AL»13l  •(l.-£KP<-0.U3*(T-TE(r)  »  » » 
ALP5*-a.  0  oi»(i.- EXP  (-o.i.j*(T-TE«  inn 
GO  TO  la 

15  ACC=ACC*OEL 

_ I  Ft.AG»I-l _ 

ALPlsALPHA ( n ’EXP (- t. 075  5* (T-rS(  1)  M 
ALP2=ALPHA(2»  ’EXPC -3.1 2*(r-T  sell  11 
ALPJ*ALPXA<TI  •E<P<-0.5»*<  T-rs/I)  )) 
ALP5«0.001*(1.-EXP«-3.12*(  r-TS(nni 
ALP10=ALP1 

GO  TO  la _ 

21  ALPlsALPHAdl 

ALP2*ALPHA(2I 
Al.PJsALPHA<  J) 

ALP5*ALPHA(5> 

GO  TO  15 

12 _ CONTINUE _ 

la  CONTINUE 

ALPSO(K<>*AI.P5 
ARG»-0£U»AU»1*:3 

IF(4RG.GT.-2tlO.)  Sl»exo(ARJI  _  _ 

Kl*KK*i 

_ <2gKK-l _ 

XSixnsXSt  KXJ  ♦A'OO'OEL 

IF(  ALPI.EO.3.1  30  TO  U  _ 

X<*(K1»*S1*(XI*(<X)  ♦A2DO*3jLI 
GO  TO  3 

X<.(K1>=XI*«>«»  ♦A70C'*3il. 

CONTINUE _ 

COMPUTE  AND  ST02E  ESTIMATED  TARGET  1/EL03ITY  AND  ACCELE^RATION 
EDM(Kl»sX3(<H -XKKII 

IF(i(2.GE.  H  ECI3M<<<«l»(E3H(<XI-E3H<K2n/'OEL 

_ X30(2>  «X3»|<1) _ _ 

IFUKK-ll  .LE.MI  30  TO  151 

00  2  IxltN02  _  _ 

J1*I*2 

AIJl.DiRA  _ 

At  Jit  Jm>PA 

CONTINUE  _ 

TmESO* (ELG-lELSt  XDEL 
C5*C0SIELG> 

Tb*-TMEBO*TAN(ELGI 
SCAL«(Ca(2l  *:?l  •*2/0EL 
C0R«C0<2»'  ,lp2*:b 

_  Aa.l)«-COP»TS _ 

A<1,2I»-;o<2»*ALP3*C3 
A«2,N01)«-C)R*A? 

AI2tNOI«A(lt2)*A2 
A(2.2l«Ta 

COMPUTE  TRANSITtON  MATRIX  EA  A<D  ITS  INTEGRA!  EA INT _ 

CALL  0SCRTIN0,3.0EL|FEtE3t'tr,5) 

CR'»«Cfl»2l*ALPA*:B 

COMPOTE  MEAN  TRACXINS  ERROR 


76 


U«ALP2*W(1I  »ALBl*i4(21+ALP^MX3(<K)*<<>IKKl) 

00  110  I<lfSn 
00  120  J>1.’40 
0(I)>0(I>«€A(I»J)>M(J) 

120  CONTlNiJE 

11 0 _ CONTINUE _ 

FClIslCB-CTlil’XKKK)  ♦CRii'Xi  <  KKI 
F»2)»XJ«KK»  •CB>:Rb*(XA(«-Xt  -X3<  KX-HI  *A2 
00  130  Isl.NO 
00  lAO  Jxtf2 

0«I)*DCI)»EAINTt I, J)*F(JI 

140 _ CJONTINUE _ 

0(I)3Q. 

130  CONTINUE 
C 

C  COMPUTE  COVARIANCE  MATRIX 

Q 

C^C  ri  .T>»«ALP5TAL’6»A0SC£3NtKKII  ♦AL*  T’ASSTEDDHt  KX»  I  > 
1  •SCAL 

CX  (2.2)*CALP5Hl«-i)*ALP%»ABS«iOH(<K-M>  ) 

1  ♦ALPT»ABS(EPOH ( X <-Nl )  1  * S CAL* A2 •• 2 
CALL  MULT  (SAINT,  COC, '10, Nl,  »l,10l 

_ C A LL  MULT(EA .p, NO. N1.P2. 13  I _ 

OO"  220  1=1. NO 
00  220  J=1,N0 
P(I.J>»P1(I,JI«'P2(I,JI 
220  CONTINUE 
150  CONTINUE 

_ OELG=ELO _ 

AZERRsHdJ/C?  * 

A2S0«S0RT(»(1,1»)/Ca 

AZTR«N<23/C3 

RETURN 

ENO 


DIMENSION  E(lItP(lltS(16tcl(ll 

00  10  1*1. L 

II»1 

00  10  K»1.L  _  _  _ 

TEMP»0, 

00  5  J«I.L1.L _ 

TEHP*TEMP*E(JI*P(III 
5  II»II*1 

XK* (K-1)*L*I 
H(KK>*TENP 
10  G«KK)«TEMP 

_ if<mr.eq.ii return _ 

00  20  1*1. L 

00  20  K*I.L  _ 

T£NP*0. 

Il’X 

00  15  J*IfLl«L 

TEMP*  TENP»G(J)*E(1H _ 

is  II*II»L 

KK*(K-1I»L*I  _  ^  , 

2J  H«<KI*TEHP 

L2*L-1 

00  30  I«1.L2 

_ L3*J*1 _ 

00  30  J«L3>L 

K1*(I-1I*L»J  _ 

K2*IJ-1I«L*I 
30  HIK1I«N(K2) 

ENO 

.^U}ROyTINE  OSCRTINOIN.A.OELtEAtEAINr.Nn 
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OIMcNSION  A(l) ,■ All) .EAin (II tC3EF<30) 

C  SETS  EAa;XB|A*0EL>t£4I'tr*INT!GRAL  EA  0  TO  OcL 

NOIHl  =  )IOIM-fl 
NN=NOIi*NOI'1 
NTH1=NT-1 

C0EF«NT>»1. _ 

do  10  I=1,NTM1 
II»NT-I 

10  COeFCIII*OEL*C3EF« II»1)/JL3AT( II 
C  IIT  MUST  a:  AT  LEAST  J 

CALL  0IA5(N3H,EAINT,A,:3E-(1I  ,;0£F(2II 

-.M  60  L^S.NT _ 

CALL  HULTIA,EAi'IT,NOtM,'IN.  EAdl 
ifil.eq.ntiso  ra  70 

6)  CALL  OIAGIMOIMdAINT.EA.l.l.COEFILII 
70  00  80  XI=1,'41|'>I3I)»1 
EA(III=EA(IT)«1.0 

80  CONTINUE _ 

ENO 

SU3R0UTINE  3TA3I NT IH,At3« 3 1» C2I 
DIMENSION  Alll.aill 
NOIHlaNOIM+l 
NN*NOIM»MOIM 

_ NMlsNDIM-1 _ 

II»1 

IFiCt.EQ.l.l)  33  TO  10 
00  5  Jst.NN.NOrM 
KsJ+NMl 
00  A  I»J,K 

_ A  A(I)aCl*9(II _ 

A(III»AIIII*C2 

5  Ilall^NQlMl 
RETURN 

10  00  7  J=1,NN,N0I1 
KaJ*NMt 

_ DO  6  IaJ,K _ 

6  AlIlaBdl 

A(  IIIaAlIIdC? 

7  II=II*N0IM1 
RETURN 

ENO 
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